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ABSTRACT: We discuss the Isgur-Wise function £(y) in the small velocity (SV) limit 
within the QCD sum rule method. The behavior of £(?/) in the SV limit is sensitive to the 
particular form of the duality relations used to decontaminate the sum rule predictions from 
the continuum contribution. Peculiarities of the duality relations in the problem at hand 
are revealed. It is shown that the proper requirements of duality and angular isotropy for 
S wave states lead to an unambiguous form of the sum rules for the Isgur-Wise function. 
We illustrate the constraints due to these requirements using a toy model of the harmonic 
oscillator. The slope parameter and the shape of £(y) are determined. 



1. INTRODUCTION. 



Strong interactions in the quark systems become simpler if one of the quarks considered 
is infinitely heavy. Then its motion is effectively classic and can be considered as given in 
the sense that the surrounding gluons and light quarks (we will refer to this component as 
to the light cloud, which is equivalent to more commonly used notion of "brown muck") 
do not affect it at all. The fact that the simplest hadronic system one can invent is built 
from an infinitely heavy quark plus the light cloud has been noted long ago . 1 The revival 
of interest to such systems we are witnessing now is due to explicit introduction of new 
symmetries 2 taking place if there are several heavy quarks, and itlq ^> Aqcd where mg 
is the heavy quark mass. The initial impetus in this direction has been given by the work 
^ where the B — > D transition has been analysed in the so called small velocity (SV) limit, 
see below. The culminating point in the recent achievements in the "heavy" quark-" light" 
quark systems is the formulation of the Heavy Quark Effective Theory 4 (for a review see 
Ref. 5), which ideally suits for describing regularities of this sector of the hadron world. 

In the semileptonic transitions of heavy mesons of different flavors the heavy quark 
symmetry is reflected in the universal Isgur-Wise (IW) function. Assume we are interested 
in the current induced transitions (it can be either vector or axial current) between 
two heavy mesons, Mj and Mf. The quark content of M^j is Qijq where Qi(Qf) is the 
initial (final) heavy quark and q is a light antiquark. For instance, we can take b as Qi and c 
as Qj. Moreover, assume for definiteness that is the ground state pseudoscalar (generic 
B). Mf can be either the ground state pseudoscalar or the ground state vector (generic D 
or D*). In general, the set of transition matrix elements of and involves a large 
number of apriori unrelated formfactors. In the limit itlq — > oo all these matrix elements 
are described by a single function which can be introduced in the following way: 

< P Qf q(v)\ vf *\ p Qd v i) >= + (l) 

Here Pq^ stands for the pseudoscalar meson with the quark content Qq, = Qfl^Qi, 
Mf(Mj) denotes the mass of the final (initial) meson, is the four- velocity, 

O = P tj/ M i,f' (2) 

and, finally, the variable y = VfV{. In the limit mg — > oo we stick to the quark and meson 
masses coincide. The peculiarity of the IW formfactor £ is that it depends only on the 
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velocity transfer ( masses are irrelevant) . If in the rest frame of the initial meson the final 
meson is also at rest then £ is trivially known 3 : 

e(y = 1) = 1. (3) 

The normalization property of the IW function in the SV limit is essentially not dynamical, 
it merely expresses the symmetry of the strong interaction. At the same time the shape 
of the function encodes dynamical properties of Quantum Chromodynamics. One can 
show that is related to the Wilson line operator for specific contours and is, thus, one 
of the most fundamental objects in QCD. Dynamical QCD-based calculations of the IW 
function are, clearly, of general interest. 

Several recent works are devoted to this issue. In Ref. 6,7,8,9 the IW function is treated 
within the QCD sum rule approach with results that seem to be quite promising. (The 
first applications of the QCD sum rules to the form factor problems are worked out in 
Ref. 11,12,13,14). Still the question is far from being closed since some aspects obviously 
require further study. One of the most evident examples is the small velocity (SV) limit 3 of 
the IW function. The slope of the function £ at y — > 1 reflects an interesting characteristics 
of the formfactor (see below). In Ref. 7 it comes out infinite, an unphysical feature of the 
corresponding analysis. A remedy suggested in Ref. 8 eliminates the spurious infinity, but 
leaves the prediction for slope parameter rather uncertain, to say nothing about an apparent 
theoretical ambiguity contained in a prescription used in Ref. 8 which calls for clarification. 

The present paper is specifically devoted to the analyses of the IW function in the SV 
limit. First, we discuss different specific aspects of the QCD sum rule calculations, elucidating 
new elements which go beyond the standard procedure 10 . The main focus is on the issue of 
how duality must be applied to estimate continuum contribution in the formfactor situation. 
In the SV limit at least one of the standard prescriptions is shown to be contradictory 
and the procedure of the continuum decontamination should be carried out in a specific 
way. Second, we comment on the additional information stemming from the Bjorken sum 
rule which relates the slope of £ in the SV limit to nondiagonal formfactors of the type 
QiQlo- ->Qf<i\i+- 

Our main results can be summarized as follows. The duality for the three-point functions 
must be understood not in a local, but in a generalized sense, as duality between integrated 
spectral densities. This is explained in detail in section 2 on the example of a toy model of the 
harmonic oscillator. The sum rules for the Isgur-Wise function and its slope parameter p 2 
are described in section 3. We argue that the symmetry requirements and the proper duality 
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lead to an unambiguous form of the sum rules. This, in turn, leads to the elimination of 
uncertainties noticed in Ref. 8,9. In section 4 we analyse the sum rules and get predictions 
for the Isgur-Wise function at small y and for p 2 . We also show the self-consistency of our 
approach using the Bjorken sum rule. 



2. CALIBRATING THE SUM RULE APPROACH 
IN THE HARMONIC OSCILLATOR. 



In this section we shall test various elements of the method in peculiar kinematics 
inherent to the problem at hand. To this end a toy model is considered where there are no 
gluons and the interaction between q and Q is pure potential. For simplicity we will choose 
the harmonic oscillator potential, 

V(f) = ^ (4) 

where r is the distance between Q and q, m is a mass of the " light quark" . 

It is convenient to work in the Breit reference frame ( Fig.l). Before the photon emission 
the heavy quark force center moves with the velocity v, after the emission its velocity is —v, 
so that the momentum transfer q=(0, q), where q = 2Mv, M — > oo is a mass of the heavy 
quark. Since our toy model is nonrelativistic (with respect to the light quark q) we have to 
assume that to Cm and, moreover, \v\ <C y/u/m. The latter constraint is imposed because 
we are interested in the SV limit. In other words the formfactor is expanded in powers of v, 
and only the terms 0(iP) and 0(-u 2 ) are kept. The first term trivially reduces to unity —an 
analog of eq. (3)— while the coefficient in front of v 2 is the central object of our analysis. 

Needless to say that in the potential model the diagonal formfactor can be found ex- 
plicitly and in the straightforward way, as well as nondiagonal transitions involved in the 
Bjorken sum rule. The nonrelativistic analogue of the IW function has the form (in the Breit 

system) 

£(^Vr. = J d\^{r)e- 2im ^{r). (5) 

Here is the ground state wave function and and the subscript "n.r." stands for "nonrela- 
tivistic" . Eq. (5) implies that for the oscillator potential (4) 



(6) 



In order to make contact with the standard parametrization accepted in the relativistic case 

^) = i-P 2 (y-i) + o((y-i) 2 ), 
y = V M ' 

it is convenient to present eq. (6) as 

f>l, = ^ (8) 

Here the fact that 

y- 1 = 2v 2 + 0(v) 4 (9) 

is explicitly taken into account. Notice that for genuinely non- relativistic systems parameter 
p 2 is large: ^ > 1. 

Along with the diagonal formfactor (5) one can consider nondiagonal transition formfac- 
tors from S to P-wave states, nonrelativistic analogues of the universal functions introduced 
in Ref. 16. In this case 

2^7-n.r = J d 3 x^j* (x)e- 2imiJ %(x) (10) 
where tpj(x) (j = 1, 2, 3) is the three-component wave function of the P-wave state: 



The explicit forms of ipQ and ipj can be found in the standard textbooks( see e.g. Ref. 17). At 
small v 2 the transition amplitude is proportional to v. In the general case there are transitions 
from the ground state ipQ to all excited P-wave states. For the harmonic oscillator, however, 
the only nonvanishing r is that due to the transition to the lowest P-wave state. 

Now, the Bjorken sum rule ^ takes the form: 

Pl.r. =rL.(v 2 = 0). (12) 

The latter equation expresses the completeness of the full set of the wave functions. From 
eq. (10) it trivially stems that 

.2 - l m 1 



T( ^ = 0) ^--^ (13) 

Comparing eqs. (8), (12) and (13) we see that the Bjorken sum rule is satisfied. Of course, 
there is no doubt that it should be satisfied. We do this simple exercise here only because we 
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want to set the framework and test the sum rule approach, to be applied in more complex 
environment of quarks and gluons. 

Let us pretend now that the wave functions are unknown, and the only quantities we are 
able to calculate reliably are the correlation functions at short (euclidean) times. Specifically, 
we consider 

S(v;n,T 2 ) = [ AK(0,ri + r 2 |f,Ti)e- Wr K(f,Ti|0,0) 

r - (14) 

This is the nonrelativistic analog of three point function. Here K(r, r|0, 0) is the amplitude 
(time-dependent Green function ) of the propagation from the point (0,0) to the point r, r 
in the Euclidean time. If t\ i — > oo only the ground state contributes and S(v, t\, t 2 ) — > 
e -E Tj^ o ^Q^|2^ nr (i? 2 ), T = r\ + T2- For the harmonic oscillator K is known exactly for all 
r (see e.g. Ref. 18), 

ts<- In n\ / mUJ \3/2 - 9 ■ ? T, ^ 2 cosh (ut) -s 
K(X, T O, 0) = ( : ) ' e 2 S inh(^r) V /_ (15) 

However we shall base our estimates on the short time expansion of K in order to imitate 
QCD calculations, and then compare the results obtained in this way with the exact ones. 

Integration over r in eq. (14) is readily carried out, and we arrive at the following v 2, 
expansion 

S(tT;Ti,T 2 ) = So + Si + ... ( 16 ) 
where Sj is the term of the i — th order in {7 2 , 

_,mw 3 , 2 1 sinh (ut{) sinh (ur 2 ) m f 2 
1 " (sinh^T)3/2 ^- 2v )■ ( 18 ) 

Notice that the zero-order term is just the same as for two-point functions (see e.g. Ref. 19), 
a straightforward consequence of the normalization theorem 3 . It is well known that the 
sum rule approach gives beautiful results for the positions of the S-wave levels in Quantum 
mechanics a refined expression of duality which works in a very transparent way in this 
case. We will concentrate on the term 0({7 2 ) determining the slope of the formfactor. 

Combining eqs. (18) and (17) in the limit T\p, ~^ 00 one would recover eq. (8) ( e.g. by 
considering the ratio Si/Sq.) We would like to examine how far one can go with the short 
r expansion. 
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Keeping only the first terms of the expansion already produces a reasonable estimate 
even disregarding the contribution due to higher states ( the latter is at the center of the 
discussion below). Indeed, in the symmetric point t\ = t<x = T/2 

Sl /S Q = {-2v 2 )--(l-±-x 2 + ...), x = coT. (19) 

CO 4 LA 

It is easy to see that the stability plateau is around x ~ 2; the corresponding estimate for 
p 2 is p 2 ~ ^(2/3) (cf. eq. (8)). At this point the power correction constitutes 1/3 of the 
leading term. Moreover, with 3 power terms included p 2 ~ (m/(2a;))0.8. Agreement with 
the exact result (8) is evident, although not very impressive, of course, which is natural, 
since in the above analysis we completely ignored contamination due to higher states. 

Let us discuss now how the higher states affect the analysis. Sq is contributed by the 
S-levels only, situated at E n = ^co + 2noo, only diagonal transitions take place, 

oo 

So = J2\MO)\ 2 e^ T , (20) 

n=0 

where 

r— r- (21) 

E n ,mw.o/ 9 2 ^ 

V CO 7T V VT 

Eq. (21) perfectly matches the "bare" spectral density, i.e. that determining So in the limit 
when the interaction is switched off, oo — > 0, 

"OO 



-L(-) 3 / 2 v^ , S bare = r o\^e- ET dE. (22) 

V2% 7T Jo 



^bare 


Indeed, the smeared sum of delta-functions, a = Y. n \Jk\J ; § i ( z ^ i ) 3 / 2 5(-E - E n ), coincides 
with eq. (22) . The standard step-like model 10 with <jq from eq. (22) and the threshold 
energy Eq = |u; gives excellent approximation for higher states 2( ^. 

Let us proceed to less studied case of S\ where, unlike Sq, we do encounter peculiarities 
of the form-factor situation. From the large r expansion of eq. (18) we learn that S\ contains 
both the diagonal and off-diagonal transitions, 

Sldiag = E(^ 2 ^)(^,,)|^(0)| 2 e-^ T 

fc=0 , (23) 

= (-2v 2 — 3 / 2 — - — - 5 / 2 -coshcuT 

v JK uo jk 2tt j K sinhuoT J 2 



and 



oo 

Sloff-diag. = E(-3 m2 ^)( r2 )fc^+l^( )^+l(°) e_i?feT ( e_2 " T2 + e_2 " Tl ) 

fc=0 (24) 

= (-2^)flffl . ' (-l/^e-^ + e-^). 
^ ^7T (sinhu/T)^ 



Here 



r 2 



| Avf (r>-* 2 ^(f), (25) 



and we have used the fact specific for the oscillator that the off-diagonal transitions are only 
from I to I ± 1 S wave state. Thus, the ground state is connected only to the first excited 
S-wave state. ipk{r) stands for the wave function of the k-th S-wave state. The set of points 
where the physical spectral density is concentrated is presented in Fig. 2. 

Now, one can easily convince oneself that the residue of the k-th state contributing to 
S\diag ( E k = (3/2V + 2wA;) is equal to 

<MS»? = -2* 2 (£)( V^^i ^(2* + 1)!}. (26) 
Furthermore, for the off-diagonal transitions, eq. (24), the residues are 
(--mV)(r 2 ) fc; , +1 ^ fc (0)^ +1 (0) = -2^ 2 (-)(— ) 3 / 2 (-l){ 3 (fcj)2 1 ^ J }■ (27) 

Notice the important minus sign appearing in passing from eq. (26) to eq. (27) . 

At this point it becomes clear that the standard picture of local duality we got used to 
in the case of the sum rules for the two-point functions 10 (see e.g. Ref. 21 ) fails in the 
formfactor situation, at least in the kinematics at hand (expansion of the formfactor near 
the SV limit). Indeed, the standard picture would imply that a few individual levels smeared 
over a natural duality domain produce (locally) the same spectral density inside this domain 
as the one obtained from free motion. We will see shortly that a conventional averaging of 
eqs. (26) and (27) leads to something quite different from the "bare" spectral density, even 
in sign. 

First we need to calculate the function S(v; t\, t^) defined in eq. (14) for the case of free 
motion {uj — 0). This is a straightforward exercise since in this case 

s ^IlW d3 ^- ¥f ^ { -^- p ^ (28) 



f is the distance between the particle and the force center; ^ and |^ are energies. 

We begin with a free particle with the momentum p + mv in the Breit frame (f? is in the 
rest frame of the force center Q, see Fig.l), and end up with the particle with the momentum 
p — mv in the Breit frame (p in the rest frame of the force center Q, after it emits a quanta 
"7" which carries away q = 2Mv, Fig.l). Since there is no interaction, nothing happens to 
the particle in the process of scattering of the force center, 

ff + mv=p — mv. (29) 

Eq. (29), obvious as it is, also formally stems from eq. (28) which can be readily transformed 
to 

Share = ^2 ^ / s dSdE'e^^ (30) 
where the integration in the right-hand side runs over the domain E depicted in Fig.3. 
In the limit t) < 1 we are interested in the bare spectral density 

1 m 

^barc(^) = (31) 

is very large and is concentrated in the narrow strip near the diagonal E = E' (Fig.3). If, 
instead of E and E', one introduces the variables 

A = E ! -E, B = -(E' + E) (32) 
2 

the boundary of E is given by the equation 

B = ^ A2 + \ m# - (33) 

Moreover, B varies in the interval [mv 2 /2, 00); for given B the maximal and minimal values 
of A are ±(B — II ^-) 1 / 2 (8mv 2 ) 1 / 2 . Thus, at v — > the width of the shaded strip in Fig. 3— 
it is proportional to v — vanishes. Naively, one would say that the bare spectral density is 
dual to diagonal transitions only (We hasten to add that the latter statement is not correct). 
The fact that <7bare > will be important in what follows. 

One more observation important for understanding how the duality can be implemented 
in the case at hand. The process with the interaction switched off which must be dual to 
that with the bound state propagation looks as follows. The propagation of the force center 
before and after scattering is like the cannon ball, (Fig.l) and we need not consider it. The 
propagation of the light "particle" is free. In the laboratory (Breit) frame we fix its energy, 
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otherwise allowing it to propagate freely in any direction it wants from the point of emission 
to the point of absorbtion. Fixing the energy of the light particle in the Breit frame means 
that we fix 

B= 1 -(E + E') = £ + ^f (34) 

where the energy in the Breit frame is denoted as £. On the other hand, allowing for all 
possible angles for the vector pq = pi + mv = p — mv implies that A = E — E' must be 
allowed to take all kinematically possible values. Thus it is clear that the domains of duality 
can be only those depicted on Fig. 4. Any other orientation of the duality domains by no 
means can be acceptable, and as a matter of fact, leads to unphysical result for the slope 
parameter as will be discussed in more detail below. 

Let us now return to the problem of duality and discuss it in more detail. Consider first 
the spectral density for the free particle. This density is concentrated in a small domain of the 
width ~ \v\ \Jm£ around the diagonal and is positive. The Taylor coefficients of expansion of 
Sb are in v 2 arise through expanding the boundary of the duality domain in v 2 and expanding 
the integration variables. The duality domain should be chosen as a certain interval in £. 
From eq. (33) we see that for fixed £ A max and A m \ n have no v 2 expansion. The only 
v 2 dependence comes from eq. (34) when one expresses B in terms of £ and v 2 and, then, 
expands in v . To illustrate this point it is instructive to consider S^are at t\ = = T/2. 
Eq. (30) implies then that 

<Wn =T2 = T/2) = —j- / d£ / dAexp -(£ + mv 2 /2)T. (35) 

8?H v Jq J Amm 



The width of the A integration is A\/2m£v . (As has been already mentioned, this width 
scales as ~ (mw) 1 / 2 t> — > at v — > 0). As a result, 



T 1 f°° 

Shared =7* = -) = -^rn 3 / 2 / d££^ 2 exp -{£ + mv 2 /2)T 
2 V27T 2 Jo 



(35) 



which coincides, of course, with the small T limit of eqs. (17) , (18) , being expanded in v 2 . 
The v 2 term in the expansion, we are most interested in, in particular, takes the form 



i r 

S 1 = -mv 2 T _ m 3 / 2 / d££ 1 ' 2 exp -£T 

2 ^ 2 t (36) 



2 1 %/o 
—mv — i=—^rn ' 
4V2tt 2 



roc 

/ d^^exp-^T. 

Jo 



The last line in this expression contains no explicit T dependence except the exponential. 
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On the other hand, consider the genuine spectral density for the oscillator. The spectral 
density for the Taylor coefficients in the expansion of S in powers of v 2n is spread over the 
length of 2nu from the diagonal (Fig. 2). The reason is that it is proportional to matrix 
elements of r 2n . For n=l we see that the spectral density is spread over the length luj as 
depicted in Fig. 2. Consequently, we see that the spectral density for S(v 2 ; t\, T2) is spread 
over the whole plane. The area in which it is concentrated has no correlation with the area 
where the free particle density is concentrated. The boundary of the area where the free 
particles spectral density is located is determined by £ and w 2 , the analogous boundary for 
oscillator is infinite (or determined by uo for Taylor coefficients). Clearly, due to this reason 
the oscillator spectral density can not be locally dual to the spectral density for free particles 
in the usual sense. We conclude that the local duality between the spectral densities for the 
oscillator and for the free particles is absent in the formfactor situation considered here. 

Although the local duality is absent, there still remains a weaker version of duality, 
namely a generalized duality. The oscillator spectral density integrated in the A-direction 
will be matched by the spectral density for the free particles integrated over the corresponding 
region, provided the weight function in the 8 direction is exponential. 

Let us explain this assertion in more detail. Start with Sq— the Taylor coefficient at zero 
order in the expansion of S in powers of v 2 . For the oscillator this density is proportional to 
\[E and is concentrated on the diagonal. For free particles this is a constant spread over a 
narrow area depicted in Fig.4. Let us now integrate the free particle spectral density between 
the boundaries in the direction orthogonal to the diagonal (here and below we shall call this 
direction the A-direction). We immediately get the spectral density ~ y/S. We conclude that 
in this case the oscillator density is dual to free particle density integrated over the permitted 
interval in A direction. Note that the latter coincides with the free particle density from the 
sum rules for the 2-point function (see eq. (23)). We thus find the origin of the duality for 
So discussed above (see eqs. (19)-(23), where we discussed the duality for Sq). Already for 
So the duality is, rigorously speaking, a generalized one. 

Consider now S\. We have already discussed above the bare spectral density. The 
integral representation for Si has the integrand ~ v and a very narrow integration domain, 
whose width in the A direction is ~ v. One can introduce an effective spectral density in 
the B (S) direction performing explicitly integration over A. Then the singularity in v is 
eliminated and the correct analytic properties in v 2 are realised in a trivial way. One can 
read off the relevant expression for the effective spectral density from eq. (36) . 

Now, proceed to the analysis of Si from the other side, with the interaction switched 
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on, oo 7^ 0. In this case our starting point is eqs. (23) , (24) , (26) , (27) . A sketch of the 
spectral density is given on Fig.2 where the k-th energy level is given by Ej~ = (3/2 + 2k)u. 
Our task is to check that the generalized duality still takes place: the appropriately defined 
sum over k matches eq. (36) . In order to check duality we must go to the limit oo / E k — > 
and check that in this case the expansion becomes that for the free particle. Here however 
we encounter a difficulty The series (23), (24) converge absolutely only when ooT is kept 
finite. Otherwise we encounter alternating series for Sdiag an d •S'non.-diag. with the terms 
increasing like /c 3 / 2 . If the contribution of the k-th diagonal transition is positive that of the 
neighboring off-diagonal transitions is negative. Generally speaking, such sum depends on 
the way of summation. The regularization of the sum is provided by the exponential factor 
exp —EfcT in eqs. (23), (24) . This factor prompts us how to carry out the summation. 

Let us consider the symmetric point t\ = = T/2, T — > (k large). It is easy to 
see that the k-th state contributing to •S'idiag nas the residue proportional to 

while that contributing to S'loff — diag nas the residue 

A; 3 / 2 + &k 1 / 2 + /3/T 1 / 2 + ... . (38) 

In terms of energies 

Suiag{E k ) ~ $>(^) 3 / 2 + a(^) 1 / 2 + ...)_L eX p (-E k T). (39) 



oo oo v w 

k 

For the nondiagonal case we get in the same way: 

5lnon.-diag.(^) ~ " £(2(% 3/2 + + ^ + -) ^ ^ E k T ) 

° * — ' 00 00 rju 

i * (40) 
(exp {—2ujt\) + exp {—2uot2}). 



'00 

Here a, (3, a, (3 are numerical constants. Notice that the weight factors in eqs. (39) , (40) 
are different. It is clear that the dangerous terms in the expansion are those proportional to 
(.E^/cu) 3 / 2 , (Efc/oo) 1 / 2 . Otherwise the series will smoothly go into (36) in the limit oo — > 0. 

Let us discuss the dangerous terms in more detail. Eqs. (23) , (24), (26), (27) can be 
combined in the following way: 

20O 7T 

k — 
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The parameter T and characteristic values of k in the sum are not independent: the choice 
of T dictates what values of k saturate the sum. The duality formally is to hold at large k 
(small T, T ~ k^ 1 + Ck~ 2 + ...). Then, one can expand the expression in braces in T, 

uj 2 T 2 

{...} = (4fc + 3) - (4k + 6)(1 - uT + _ + ...) (42) 
= -3 + 4kuT + ... . 

Notice that in 0(T®) the terms ~ k cancel. This means that under the appropriate summa- 
tion the terms ~ (E^/uj) 3 / 2 in eqs. (39) and (40) annihilate each other. 

The next step is to inspect what happens at the next-to-leading order. Since the indi- 
vidual contributions grow like k 3 / 2 it is quite obvious that the sum is saturated by terms 
with k = 4^ (one can see this from maximizing k 3 / 2 exp (— 2kuT). As a result, —3 in eq. 
(42) is cancelled by AkcuT. In terms of eqs. (39) and (40) this fact implies cancellation of 
0((E k /u) l l 2 ) parts. 

As a result we are left with terms O^/cu)" 1 / 2 ). To see that the coefficient matches 
that in eq. (36) one must keep the term T 2 in eq. (42) as well to use a relation between T 
and k, T = ^ + Ck~ 2 , truncated at 0(k~ 2 ) level. 

Notice that our way of summation of diagonal and nondiagonal parts precisely corre- 
sponds to the integration in the A direction. We conclude that the spectral density integrated 
in the A direction for the oscillator is dual to the spectral density for free particles, also in- 
tegrated in the A direction. Notice that duality is valid not for local densities cr(E, E ), but 
only for integrated densities in the A direction er(£). It is also obvious, from the discussion 
above, that the integration in the A direction is the only way to ensure duality. Otherwise 
we encounter dangerous terms ~ 1/ou and/or can even get terms with the wrong analytic 
properties ( not expandable in v 2 )7 

These considerations directly affect the possible choice of the continuum model. It is 
clear that the considerations above lead to the choice of the continuum model depicted on 
Fig. 4 as the only possibility . 

We discussed above an analogue of the double Borel transformed sum rules. During 
the investigation of the 3-point functions some authors considered also the so called sum 
rules in the external field 22 >23, 13,24^ ^at correspond to single Borel transformation. A 
non-relativistic analogue of the single- Borel transformed QCD sum rules is the function 

S sm = J2(S)nn\M0)\ 2 exp -E n T - Sm ^ ( "^ } ~ ^"^V n Wf(0)- ( 43 ) 



n n ,l 
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Here 

S nl = J d 3 f^ n (r) exp (-2imvr)^ (f) . (44) 

It is clear that this function differs from the function S given by eq. (14) (taken at t\ = tq) 
by different weight for the nondiagonal transitions. The single Borel sum rules (we shall use 
this name instead of more lengthy single time Borel transformed sum rules) are widely used, 

OO O^i 1 'X 0/1 

e.g. in the study of the pion and proton formfactors ^ z °i 1J ^. This is the reason why we 
consider them here. It is easy to see that in the limit T — > oo 

S™{v\T) - (\m\ 2 u,.(v 2 ) + ^pV^T (45) 

Here const = — ipo(0)ipl(0)Soi + ... measures the contribution of the nondiagonal transitions. 
Asymptotically at T — > oo 

uTS sm (v 2 , T) = (a x + b)e- E » T (46) 

where x = uoT and a = |^(0)| 2 £ n . r .(v 2 ). From eqs. (23) , (24) it is easy to find exactly the 
corresponding coefficient for the terms ~ v 2 in S sin : 

„^m.mw.o/ 9 1 1 sinhcjT. 

The r.h.s. of eq. (47) multiplied by (sinh (wT)) 3 / 2 wT is well approximated by a straight line 
(46) for x = ujT for sufficiently large x. Indeed, it has the form 

Sf^T _ ( .ml 



-2^) (xcoth(x) - 1). (48) 



5 2 
Taking several first terms in the expansion around point x = Owe have: 

^= (-2^)^/6 + ...). (49) 

The r.h.s. of eq. (49) is well approximated by a straight line for x > 2. We get p 2 ~ 
Taking 3 power correction terms we get p 2 ~ ^0.8. This can be compared with eq. (19) 
and the discussion that followed. We see that in the case of the harmonic oscillator the single 
Borel sum rule works as good as double Borel sum rule. 

We note however, that in the case of the electromagnetic formfactor of pion the single 
Borel sum rules give overestimated value of the formfactor and underestimated value of the 
slope parameter, although within the error bars (see Ref. 22,23 ). 

14 



Consider now the duality for the single Borel sum rule, i.e. function S sm . Sjdiag i s 
the same as in eq. (23) . The off-diagonal part is almost the same as in eq. (24) , with 
n = r 2 = T/2. The only difference is that e~ 2wTl + e~ 2uJT2 = 2e~ wT in this equation is 
substituted by (1 - e" 2wT )/(cjT) in eq. (47) . These two expressions have identical first two 
terms in the loT expansion. In accordance with the discussion above this means that the 
dangerous terms in S sm are precisely the same as the dangerous terms in the double Borel 
sum rule considered above. The duality considerations are the same for both sum rules. 
When we write this sum rule in terms of the spectral density cr(S) we do it after implicit 
integration in A direction. 

Let us now summarize lessons that follow from the study of duality for three point func- 
tions for harmonic oscillator. Note that these lessons refer not only to the case at hand 
discussed in the remainder of this article , but also to other cases of the calculation of 
three-point functions in QCD in similar kinematics. First, considerations of duality and 
appropriate symmetries impose strict constraints on the continuum model which must cor- 
respond to the region of Fig. 4 . Second, duality is not local. The duality in three- point 
functions is between integrated in the A direction spectral densities. 



3. SUM RULES FOR THE ISGUR-WISE FUNCTION. 



In this section we shall write the sum rules for the IW function and extend the analysis 
of the preceding section to the heavy quark case. 

We begin from the double Borel sum rules for the IW function in the symmetric point 
°\ = °2 = 2cr, where a denotes the Borel parameter. This sum rule was already written in 
Ref. 6,7,8,9. We refer the reader to these articles as well as to Ref. 1 for the general discussion 
of Borel transformed nonrelativistic sum rules in heavy quark systems and the derivation 
of the sum rule for the IW function. We have rechecked that the coefficients in front of 
condensate terms are correct. The bare spectral density can be obtained from calculations 
in Ref. 25 which were carried out for finite heavy quark masses. Using explicit expressions 
in that work it is not difficult to take the limit mq — > oo. 
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The sum rule for the IW function has the form: 



3 I I A A I \ ( S l + S 2) 

1 ' ds\ds2\S\ + S2) exp 



&7r 2 (y + l)y/y 2 -lJ h " ' 2cr 

6(2y+l) F 2 S(y)exp-E R /a 
a/ + 263^2 " 4 

Here a, 6 are the vacuum condensates averages: 

a =< 0(0)0(0) > = -(0.24) 3 GeV 3 ; 

b =< ^(0)^^(7^-0(0) >o= ^o a ; m l ~ °- 8 GeyS 

y is the recoil parameter, a is the Borel parameter. 



(50) 



(51) 



Our parameters are connected with the analogous parameters used in Ref. 8 in the 
following way. The Borel parameter T used in Ref. 8 is T = 2a. The spectral density 
variables used in our eq. (50) are connected with the dimensionless spectral variables used 
in eq. (3.14) in Ref. 8 as s\ = zcr, S2 = z'a. The continuum threshold in that paper oo c is 
numerically twice as big as the continuum threshold oo c used in our paper, since the author 
uses the spectral density variables 00,00' equal to 2s\, 2s2 respectively. 

Recall that in the nonrelativistic case we make the Borel transform over the light quark 
virtual energy and a has the dimension of mass, not of mass square as in the relativistic case. 
This makes the sum rules very similar to the sum rules for the harmonic oscillator considered 
in the previous section. The area of the integration of the spectral density is denoted by £ 
and is depicted in Fig. 5. Finally, Er is the mass difference between the mass of the heavy 
quark and that of the heavy meson (both tend to infinity in our limit). The constant F is 
the heavy meson leptonic decay constant: 

iF 

<0\q lfll5 h Q \P >= -^=Vf, (52) 

The field Hq are heavy quarks fields in Heavy Quark Effective Theory 4 . The constant F is 
related to more familiar fpbjF = fpyfrriQ modulo logs of the type discussed in Ref. 26 . 

Although the sum rule is the same, our analysis will be quite different from the analysis 
in Ref. 7,8,9, however. In particular we shall argue that the uncertainty in p 2 j s quite small, 
not gigantic as in the latter papers. Indeed, the main unsolved problem in the analysis 
of the sum rule (50) is the choice of continuum. We saw in the previous section that the 
considerations due to duality and relevant symmetry severely constrain the possible choice 
of continuum model in the sum rules. The analysis in the previous section can be exactly 
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repeated for the nonrelativistic sum rules in QCD. As a result the only possible choice of the 
integration area is that depicted in Fig. 5. For every fixed value of (si + S2)/2 we have to 
integrate over all possible values of A = (s\ — 82)- Any other prescription for integration of 
spectral density will be incompatible with the general requirements of duality and angular 
isotropy (S wave states). The variable (si + S2}/2 plays the same role as £ in the previous 
section. 

We have argued in a very similar kinematics that duality in the three point functions 
must be understood in the generalized sense. Namely, the spectral density for the interacting 
case integrated in the A direction matches that for the free particle only after integration 
in the A direction. The same kind of generalized duality must also hold for the QCD sum 
rules. 

The generalized duality means that we have to write the sum rule for the spectral density 
already integrated in the A direction (recall that this is the direction orthogonal to the main 
diagonal). This effective spectral density will be dependent only on the variable (s\ + S2)/2. 
Only after this integration we can speak about duality and make a choice of continuum 
model. We now have the only one possible choice of the continuum model. The continuum 
is approximated by the spectral density due to the bare loop, integrated over the variable 
A — {s\ — S2). This is exactly the same continuum model as the model that was used in 
the two-point sum rules. Numerically, of course, it does not matter whether we assume 
local duality and take the bare loop contribution integrated over a two-dimensional area 
depicted on Fig. 5 as the continuum model or if we approximate the A integrated continuum 
contribution by the bare spectral density, integrated in the A direction, starting from some 
threshold. However our experience with the harmonic oscillator shows, that only the second 
prescription is physically correct. 

We can now write a sum rule for the three point function based only on the generalized 
duality discussed above. It has the form: 



1^ s 2 e- s '°ds --+ K2y + l) = ^(^W-^A 7 
Jo 4 192a 2 4 



( 2/ + i)2 7r 2y — - 4 1 192a 2 4 • ( 53 ) 

Here s — (s± + S2)/2. This sum rule coincides with one of many possible forms of the sum 
rules (due to different choices of the integration domain E in eq. (50)) considered in Ref. 8. 
Now we know that eq. (53) is the only possible form of the sum rule. We also write the sum 
rule for the slope parameter p 2 directly: 

Jo 6 96a 2 j F 2 /4 9 ' ^ ] 



4tt 2 
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Consider now the single Borel sum rule. This sum rule was found to work as well as the 
double Borel sum rule in the case of the harmonic oscillator. We test this sum rule in the 
heavy quark case. 

Recall the basic idea behind this approach. We write the dispersion relation for the 
polarization operator 

mM=r lma ^: S2 f F j s (55) 

JO (si + Ei){s 2 + E 2 ) 

where 

(vf + v$)n{E 1 ,E 2 ,y) = J d A xd A yex P i(k lX + kz) < j 5 (x)j^z)j+(0) > . (56) 

In eq. (56) k\ is the momentum of the initial heavy meson, k is the transferred momentum, 
k 2 — k\ + k is the momentum of the final heavy meson; k^ = (mg + Ei)vi^ and k 2fl = 
(rriQ + E 2 )v 2fJj . The current j'5 = Qj^q is the quark current that creates heavy mesons, and 
j > — Ql{iQ is a vector current. Q and q are the heavy and light quark fields respectively. 
In the double Borel transformed sum rules we make the double Borel transformation of the 
sum rule (55) with respect to the the variables E\ and E 2 independently and saturate the 
r.h.s. of eq. (55) by the resonance contributions. In the single Borel transformed sum rules 
we put E = E\ = E 2 from the very beginning in eq. (55) and then carry out the Borel 
transformation over the variable E. As it was explained in the previous section in the example 
of the harmonic oscillator the difference between the double and single Borel transformed 
sum rules is that we take the contribution of nondiagonal transitions with different weights. 
The spectral density for the bare loop now depends only on one variable E. As it was 
explained above when we use such spectral density we do it after an implicit integration 
in the A direction. Thus, the two sum rules are very similar in spirit. The single Borel 
transformed sum rule has the form: 

U(x, y, u) c ) = + a(y)a. (57) 

Here a(y) measures the contribution of nondiagonal transitions. One can easily calculate 
II in this case taking into account the same operators as we took into account writing the 
double Borel sum rule. The result is 

i^T (1 + {V + ^ f (58) 

- «/4 + = ^(t(y) + a(y)a) exp (-E R /a) 
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4. ANALYSIS OF THE SUM RULES. 



and 



Let us now analyse our sum rules. We start from the sum rule (53). In order to minimize 
the effects related to the uncertainty in F, we deal with the sum rule (53) divided by 
the sum rule for F,E R (Ref. 1,7,8,9 ) : 

m = I -^ 2 H = G(a,y,u, c ). (59) 
1m{v, Uc) 

Here 

3 [ u ° o ,/rr , a b(2y+ 1) 

'^.-cH^/^V^-^ + ^-j. (61) 

The function Ij^(a,u} c ) is the theoretical part of the sum rule for the mass and leptonic 
decay constant of heavy meson. 

F 2 

I M (a, u c ) = — exp (—Eft/ a) (62) 

The continuum threshold is chosen to be the same as for the two-point sum rules. We did 
not find any significant dependence of the continuum threshold on the recoil y. We depict 
the behaviour of £(y) for different values of the continuum threshold in Fig.6. We find only 
a small uncertainty ~ 10% due to the choice of the continuum threshold. Let us also remind 
once again that we assume that the generalized duality picture is valid. Hence our continuum 
model looks like the continuum model for the sum rules for the two point function and has 
a single threshold u c . We vary this threshold to obtain the best stability of the sum rules. 

Let us now proceed to the sum rules for p 2 , eq. (54). The dependence of the r.h.s. of 
this sum rule on a is depicted in Fig. 7. We immediately get: 

p 2 = 0.7±0.1. (63) 

Consider now the single Borel transformed sum rules. We rewrite this sum rule in the 
form: 

<y)a + «») = (64) 

i M (cT, UJ C ) 

Here the function / is the l.h.s. of the sum rule (58). The function a(y) is the x— independent 
function that measures the contribution of the nondiagonal transitions. The l.h.s. of the 
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sum rule (64) is well approximated by a straight line for a > 0.8 GeV. The corresponding 
behaviour of the Isgur-Wise function is depicted in Fig.8. Notice that the sum rule (60) 
gives slightly smaller (~ 10% at y — 2) values of than the sum rule (54) . Nevertheless, 
the difference between the predictions of these two sum rules is inside error bars for each of 
them. 

We can now differentiate the sum rule (60) with respect to y at y — 1. Then we get the 
sum rule for p 2 . This sum rule gives the value 

p 2 = 0.55 ±0.1 (65) 

Thus we get a slightly smaller value of p 2 using the sum rule (60). The difference is, however, 
within the error bars. 

It is tempting to speculate that the true value of p 2 lies in the intersection of the predic- 
tions of the single and double Borel sum rules. After all, these two sum rules worked equally 
well in the case of the harmonic oscillator. This gives p 2 ~ 0.6 — 0.65. 

For more careful estimate, we take the average of the predictions of these two some rules. 
We have 

p 2 = 0.65 ±0.15 (66) 
as the average from the predictions of these two sum rules. 

Let us now turn to the check of self-consistency of our approach. The nontrivial check 
is provided by comparison of our results with the predictions of the Bjorken sum rule 15 . 

The Bjorken sum rule 15 has the form 

-i m max Pmax 

^ = 4+ £ \t™ /2 (1)? + 2J2K /2 (1)\ 2 + .... (67) 
m=l to=1 

Here "/"-^(l), 73/2(1) are nondiagonal transition formfactors from the multiplet of pseu- 
doscalar and vector mesons into multiplets with sj l = l/2 + and 3/2 + respectively taken at 
zero recoil. Here si and 717 are the spin and parity of the light quark in the meson. These 
formfactors are defined in the following way: 



< A\A fi \P> 
m 



y=l = - 4t -i/2K - v 2)» (68) 



m 



y=l = ~~yf ir 3/2 W7 e a ( Vl + v ^ ^ ~ v ^ ( 69 ) 
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Here A is the scalar heavy meson in 1/2+ multiplet, and B denotes the state with J p = 1+ 
in 3/2 + multiplet of heavy mesons. We use the notations of Ref. 14,16. For general values 
of recoil y the formfactors r-^, r 3 / 2 are new universal functions of the recoil analogous to 
the Isgur-Wise function. All other formfactors for transitions between l/2~ and 1/2+ , 3/2 + 
multiplets can be expressed through r^ 2 an d r 3/2- 16 

It is possible to analyse the formfactor 71/2(1) using the QCD sum rules in which we do 
double Borel transformation. Taking into account only the contribution of the unit operator 
and the quark condensate we obtain the following sum rule for T\j2- 

lim — 3 / (si - si) exp {-s\jo\ - s 2 / 1 ai)ds\ds<i- 

a/2 = 8r 1/2 ^P exp {-E s /a x - E P /a 2 ) 

A rough estimate using this sum rule gives 

r 1/2 (l) ~ 0.35 ±0.2 (71) 

Note very high uncertainty in this estimate. First, we took into account only two operators of 
dimensions zero (unit operator) and 3 (light quark condensate). Second, and more important, 
we did not take into account continuum contribution. In eq. (70) A is a semileptonic decay 
constant of the corresponding resonance: 

< 0\j\S(P) >= \ s{P) <i>s(P) (72) 

where j is the quark current that creates a meson , Qq for scalar meson S and Q75<? f° r 
pseudoscalar meson P, is the corresponding wave function. 

It is interesting to note that the same ratio Ti^/^l) can be obtained by combining 
nonrelativistic quark model and predictions obtained using the nonrelativistic sum rules for 
the two-point functions. 1 Indeed, in the nonrelativistic quark model 16 

ri/2(l) ~ J dhtPl^kfaoik) (73) 

Here <pn is the radial wave function of the scalar resonance, and 0oo - °f pseudoscalar meson. 
For the rough estimate we have 



m T mW \4>n(0)\ Mp „, 

T SI % M <74) 
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The latter equation follows from 1 

A 2 = 12|0(O)| 2 



(75) 



The determination of T3/2 is much more cumbersome since the corresponding currents 
contain derivatives. However, following the nonrelativistic quark model we can accept for 
orientation that T3/2 ~ r i/2- Note also that the contribution of the radially excited states 
to Bjorken sum rule is strongly suppressed (it is zero in nonrelativistic quark model). We 
now take for r-jy 2 the value (see eqs. (71), (74)) T-jy 2 (l) ~ r 3/2(l) ~ 0.35. Then we get for 
p an estimate: 

p 2 ~ i + 3r 1/2 (l) 2 ~ 0.25 + 3(0.35 ± 0.2) 2 (76) 
We finally conclude that Bjorken sum rule implies 

0.35 < p 2 < 1.15 (77) 

This result is in agreement with the sum rule prediction (63) . 



5. CONCLUSION 



We have used the QCD sum rule method to find the shape of the Isgur-Wise function 
and its slope parameter p 2 . We have discussed the nature of the uncertainties pointed out in 
Ref. 8,9 and shown that they can be eliminated if we implement the general requirements of 
duality in the appropriate way. We used a toy model of a harmonic oscillator to discuss in 
detail the consequences of the demand of duality and relevant symmetries for the sum rules 
for formfactors. 

In particular, we have shown on the example of harmonic oscillator that there is no local 
duality in the sum rules for the three point functions. Duality must be understood in a 
generalized sense, between integrated quantities. The proper duality severely constrains the 
form of the continuum model. 

For the case of heavy mesons the requirements of duality lead to the double Borel sum 
rules (53), (54). The corresponding behaviour of the Isgur-Wise function as a function of 
recoil is depicted in Fig. 6. For the slope p 2 we get 

p 2 = +0.7±0.1 (78) 
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The single Borel transformed sum rules lead to the shape of Isgur -Wise function depicted 
on Fig. 8. They lead to the slightly smaller (within 10%) value of The difference 

between the predictions of these two sum rules is within the error bars however. For the 
slope parameter we obtain 

p 2 = 0.55 ±0.1 (79) 

For p 2 we can argue that the values (79) and (78) are upper and low bounds on the value of 
the slope parameter. Taking the average of these two estimates , we have 

p 2 = 0.65 ±0.15 (80) 

Note that our predictions for the slope parameter lead to the smaller values than those 
obtained by analysing the phenomenological data: p ~ 1.4±0.4 in Ref. 27 and p ~ 1.2±0.4 
in Ref. 28. 

On the other hand, our results for the slope parameter are very close to the results of 
Ref. 14. The model developed in the latter paper gives the value p 2 ~ 0.7. 

After this work was completed we have learnt about the results for p 2 obtained from 
the so called optical sum rules 29 extending the Bjorken approach. 15 The analysis in Ref. 29 
leads to the inequality p 2 < 0.75. This is in a very good agreement with our estimate for p 2 
given by eqs. (78) , (79), (80). 

We are grateful for to A. Vainshtein and I. Kogan for useful discussions. 
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Figure Captions 

Fig.l: Heavy quark in Breit reference frame a) before and b) after the emission . 

Fig. 2: The spectral density for p 2 r for the oscillator. The difference between • and 
is the sign of the residue. 

Fig. 3: Boundaries of the domain where the spectral density Ob are ^ 0. 

Fig. 4: The only possible geometry of the duality domains. 

Fig. 5: The duality domain for the sum rule for Isgur-Wise function. 

Fig. 6: The Isgur-Wise function i(y). 

Fig. 7 The dependence of the l.h.s. of the sum rule for p 2 on the Borel parameter x. 
Fig. 8 The Isgur-Wise function from the single Borel transformed sum rule. 
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